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in  deployed  US  military  service  members: 
prospective  cohort  study 

George  A.  Bonanno,  Anthony  D.  Mancini,  Jaime  L.  Horton,  Teresa  M.  Powell,  Cynthia  A.  LeardMann, 
Edward  J.  Boyko,  Timothy  S.  Wells,  Tomoko  l.  Hooper,  Gary  D.  Gackstetter  and  Tyler  C.  Smith, 
for  the  Millennium  Cohort  Study  Team 


Background 

Most  previous  attempts  to  determine  the  psychological  cost 
of  military  deployment  have  been  limited  by  reliance  on 
convenience  samples,  lack  of  pre-deployment  data  or 
confidentiality  and  cross-sectional  designs. 

Aims 

This  study  addressed  these  limitations  using  a  population- 
based,  prospective  cohort  of  US  military  personnel  deployed 
in  support  of  the  operations  in  Iraq  and  Afghanistan. 

Method 

The  sample  consisted  of  US  military  service  members  in  all 
branches  including  active  duty,  reserve  and  national  guard 
who  deployed  once  (n  =  3393)  or  multiple  times  (n  =  4394). 
Self-reported  symptoms  of  post-traumatic  stress  were 
obtained  prior  to  deployment  and  at  two  follow-ups  spaced 
3  years  apart.  Data  were  examined  for  longitudinal 
trajectories  using  latent  growth  mixture  modelling. 

Results 

Each  analysis  revealed  remarkably  similar  post-traumatic 


stress  trajectories  across  time.  The  most  common  pattern 
was  low-stable  post-traumatic  stress  or  resilience  (83.1% 
single  deployers,  84.9%  multiple  deployers),  moderate¬ 
improving  (8.0%,  8.5%),  then  worsening-chronic  post- 
traumatic  stress  (6.7%,  4.5%),  high-stable  (2.2%  single 
deployers  only)  and  high-improving  (2.2%  multiple  deployers 
only).  Covariates  associated  with  each  trajectory  were 
identified. 

Conclusions 

The  final  models  exhibited  similar  types  of  trajectories  for 
single  and  multiple  deployers;  most  notably,  the  stable 
trajectory  of  low  post-traumatic  stress  pre-  to  post¬ 
deployment,  or  resilience,  was  exceptionally  high.  Several 
factors  predicting  trajectories  were  identified,  which  we  hope 
will  assist  in  future  research  aimed  at  decreasing  the  risk  of 
post-traumatic  stress  disorder  among  deployers. 
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The  operations  of  the  US  military  since  2001  have  raised 
important  concerns  about  the  psychological  costs  of  intense  and 
repeated  deployments. 1-3  Military  personnel  deployed  for  combat 
operations  are  frequently  exposed  to  traumatic  events  and, 
predictably,  those  returning  from  the  operations  in  Iraq  and 
Afghanistan  have  evidenced  increased  rates  of  post-traumatic 
stress  disorder  (PTSD).1’2  Unfortunately,  it  is  difficult  to 
accurately  assess  the  psychological  cost  of  deployment  in  these 
samples.4  Although  previous  studies  have  made  important 
advances  in  our  understanding  of  deployment-related  trauma 
reactions,  with  few  exceptions5,6  little  prospective  research  among 
population-based  samples  exists  to  explain  how  trauma  symptoms 
vary  across  time.  These  earlier  studies  have  typically  relied  on 
convenience  samples,  which  potentially  introduce  sampling  bias 
and  limit  generalisability;  use  of  assessment  data  that  may  be 
vulnerable  to  reporting  bias;  and  cross-sectional  research  designs 
that  cannot  distinguish  trauma  reactions  from  pre-existing 
symptoms. 

The  current  investigation  attempted  to  fill  this  knowledge  gap 
by  using  data  from  a  large  population-based  sample  of  military 
personnel,  the  Millennium  Cohort  Study,  followed  prospectively 
over  nearly  a  decade.7  A  previous  study  using  data  from  this 
cohort  identified  new-onset  PTSD  symptoms  in  7.6-8.7%  of 
deployers  who  reported  combat  experiences,  compared  with 
1. 4-2.1%  of  deployers  who  did  not  report  combat  experience.2 
What  is  not  yet  known,  however,  is  how  post-traumatic  stress 
varies  from  pre-  to  post-deployment,  the  extent  that  different 
trajectories  of  post-traumatic  stress  might  capture  that  variation 


and  the  predictors  that  might  best  inform  those  trajectories.  To 
address  these  questions,  we  examined  pre-post  trajectories  of 
post-traumatic  stress  from  participants  who  had  deployed  at  least 
once  in  support  of  these  operations. 

To  identify  the  trajectories,  we  employed  latent  growth 
mixture  modelling  (LGMM).  Traditional  statistical  approaches 
used  to  examine  mental  health  over  time  assume  a  single 
homogeneous  distribution  of  change.8  Recent  advances  suggest 
that  responses  to  potential  trauma  are  heterogeneous  and 
explained  by  a  finite  set  of  longitudinal  trajectories,9  including 
chronically  elevated  post-traumatic  stress  following  the  event, 
continuous  symptom  elevations  that  preceded  the  event  and 
elevated  pre-deployment  (baseline)  symptoms  followed  by  steady 
improvement.  This  previous  research  has  also  repeatedly 
demonstrated  resilient  outcomes  characterised  by  a  stable 
trajectory  of  healthy  adjustment  (for  example  little  or  no 
symptoms  pre-  to  post-event).10  Latent  growth  mixture  modelling 
is  well-suited  to  capture  these  trajectories  because  it  explicitly 
assumes  outcome  heterogeneity.11  Of  note,  LGMM  has  been  used 
effectively  to  identify  prototypical  trajectories  following  different 
types  of  potential  trauma,  including  traumatic  injury,12  natural 
disaster,13  disease  epidemic,14  bereavement15  and  life-threatening 
medical  procedures.16  A  recent  study  mapped  trajectories  of  PTSD 
symptoms  prospectively  in  a  sample  of  US  soldiers  deployed  on  a 
6-month  NATO-led  peacekeeping  mission.17  Four  distinct 
trajectories  were  found  in  this  sample,  with  84%  categorised  in 
a  resilient  pattern  of  low-stable  post-traumatic  stress  from  pre- 
to  post-deployment.  Surprisingly  absent  was  a  chronic  PTSD 
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pattern,  which  was  attributed  in  part  to  the  fact  that  peacekeepers 
experience  less  exposure  to  war-related  trauma  than  is  normally 
the  case  for  combat  veterans. 

Pre-deployment  and  follow-up  data  from  the  Millennium 
Cohort  Study,  the  largest  longitudinal  military  study,  were  used 
to  investigate  the  trajectories  of  post-traumatic  stress  symptoms 
among  service  members  prior  to  and  after  deployment. 


Population  and  data  sources 

The  Millennium  Cohort  Study  began  collection  of  self-reported 
health  data  in  2001,  just  prior  to  the  start  of  military  operations 
in  Iraq  and  Afghanistan.  Invited  participants  to  the  first  panel 
were  randomly  selected  from  US  military  personnel  serving  in 
October  2000,  oversampling  for  women;  those  previously 
deployed  to  Bosnia,  Kosovo  or  Southwest  Asia;  and  reserve  and 
national  guard  members.  Participants  are  surveyed  every  3  years 
throughout  a  21 -year  planned  follow-up  period.  Detailed 
descriptions  of  methodology  for  the  Millennium  Cohort  Study 
are  available  elsewhere.7'18 

Of  the  77  047  consenting  participants  who  enrolled  in  the 
study’s  first  panel  (36%  baseline  response  rate),  55  021  (71%) 
completed  the  first  follow-up  (2004-2006)  and  54  790  (71%) 
completed  the  second  follow-up  questionnaire  (2007-2008).  The 
Defense  Manpower  Data  Center  supplied  demographic  and 
deployment-related  data  obtained  from  electronic  personnel  files. 

This  study  examined  the  first  panel  of  participants  who 
completed  all  three  questionnaires  as  two  separate  subpopulations: 
(a)  those  who  had  their  first  complete  deployment  in  support  of 
the  operations  in  Iraq  and  Afghanistan  between  baseline  and  the 
first  follow-up  and  did  not  deploy  again  before  April  2009  (single 
deployers);  and  (b)  those  who  had  their  first  deployment  in 
support  of  the  operations  in  Iraq  and  Afghanistan  between 
baseline  and  the  first  follow-up  and  deployed  at  least  once  more 
in  support  of  the  same  operations  thereafter  (multiple  deployers). 

This  research  has  been  conducted  in  compliance  with  all 
applicable  federal  regulations  governing  the  protection  of  human 
subjects  in  research  (Protocol  NHRC.2000.0007). 

Primary  outcome 

Post-traumatic  stress  was  assessed  at  all  three  time  points  using 
the  PTSD  Checklist  -  Civilian  Version  (PCL-C).19  The  PCL-C  is 
a  17-item  self-reported  measure  of  PTSD  symptoms  requiring 
participants  to  rate  the  severity  of  each  symptom  during  the 
previous  30  days  on  a  Likert  scale.  Post-traumatic  stress  scores 
were  examined  continuously,  ranging  from  17  to  85  points  (higher 
scores  indicating  more  symptoms  or  greater  severity). 

Covariates 

The  following  covariates  were  assessed  pre-deployment  (baseline). 
Using  the  Short  Form  36-Item  Health  Survey  for  Veterans  (SF-36V), 
the  Physical  Component  Summary  (PCS)  score  was  used  as  a 
continuous  measure  to  assess  physical  health.20  Higher  scores 
reflected  more  favourable  health  status.  History  of  life  stress  was 
assessed  as  a  continuous  measure  by  applying  the  points  and 
scoring  mechanisms  of  the  Holmes  and  Rahe  Social  Readjustment 
Rating  Scale.21  Each  affirmative  response  to  a  stressful  event,  such 
as  divorce,  major  financial  problems  and  violent  assault,  was 
weighted21  and  summed,  so  higher  scores  reflected  higher  stress 
levels.  Two  binary  variables  were  used  to  assess  alcohol 
consumption.  Men  who  reported  drinking  more  than  14  alcoholic 
drinks  per  week  and  women  more  than  7  drinks  per  week  were 


considered  to  be  heavy  drinkers.  Men  consuming  5  or  more  drinks 
per  occasion  or  day  and  women  who  reported  4  or  more  were 
considered  binge  drinkers.22  Based  on  questions  that  assessed 
lifetime  smoking  habits,  participants  were  classified  as  never 
smokers,  past  smokers  or  current  smokers.23 

Combat  experience  was  assessed  as  a  score  from  0  to  5  based 
on  the  number  of  affirmative  responses  to  personal  exposure  in 
the  3  years  before  response  to  the  follow-up  questionnaire  to: 
witnessing  a  person’s  death  due  to  war,  disaster  or  tragic  event; 
instances  of  physical  abuse  (torture,  beating,  rape);  dead  and/or 
decomposing  bodies;  maimed  soldiers  or  civilians;  and  prisoners 
of  war  or  refugees.3  Injury  status  was  derived  from  the  Joint 
Theater  Trauma  Registry,  a  database  established  in  2002  and 
maintained  by  the  US  Army  Institute  of  Surgical  Research.  It 
contains  medical  information  on  military  personnel  who  have 
sustained  either  non-combat  or  combat-related  injuries  during 
deployment.  For  single  deployers,  participants  were  categorised 
as  injured  if  an  injury  was  recorded  in  the  Joint  Theater  Trauma 
Registry  between  the  baseline  and  first  follow-up.  All  other 
participants  were  categorised  as  injured  if  an  injury  was  recorded 
between  their  baseline  and  second  follow-up  questionnaire. 

Statistical  analyses 

Separate  LGMM  for  single  and  multiple  deployers,  using  the 
Mplus  statistical  modelling  programme,  was  conducted.11  Latent 
growth  mixture  modelling  began  with  a  single-class  growth  model 
and  then  added  classes  to  examine  improvements  in  model  fit.24’25 
The  final  model  solution  was  determined  by  a  combination  of  fit 
statistics,  conceptual  rationale  and  interpretability.1 1,24,25 

We  examined  each  data-set  separately  in  three  steps.  First,  we 
identified  a  single-class  growth  model  without  covariates.  Second, 
we  compared  one-  to  five-class  unconditional  LGM  models  (no 
covariates)  by  assessing  relative  fit  with  conventional  indices 
(Table  1),  which  penalise  more  complex  models  when  they  fail 
to  provide  a  better  fit  to  the  data.24-26  To  ensure  that  the  restricted 
populations,  which  excluded  participants  with  missing  data,  were 
representative  of  the  larger,  unrestricted  populations,  we 
compared  the  unconditional  models  of  the  restricted  populations 
with  the  unrestricted  populations.  Third,  consistent  with 
recommendations  for  correct  model  specification,  the  LGMM 
was  extended  to  include  predictors  of  class  membership  in  a 
conditional  model.  This  step  is  increasingly  recommended,24  since 
the  absence  of  significant  relationships  between  theoretically 
relevant  covariates  and  class  membership  can  indicate  an  incorrect 
local  solution.  As  recommended,11  we  examined  covariates  based 
on  their  conceptual  relevance  and  performance  in  preliminary 
univariate  testing.  Finally,  multinomial  logistic  regression  analyses 
were  nested  within  the  LGMM  to  examine  predictors  of  class 
membership. 

Results 

Of  the  77  047  participants,  46  438  completed  both  follow-up 
questionnaires.  Among  those  participants,  7272  had  missing 
outcome  information  and  2946  had  incomplete  covariate  data 
and  hence  were  excluded,  leaving  36  220  participants.  Of  these 
eligible  participants,  there  were  7787  participants  who  had  their 
first  complete  deployment  between  their  baseline  and  first 
follow-up  questionnaire.  For  3393  of  these  participants,  this  was 
their  only  deployment  until  April  2009  (i.e.  single  deployers). 
The  remaining  4394  participants  deployed  one  or  more  times 
thereafter  (i.e.  multiple  deployers). 

Characteristics  of  single  deployers,  multiple  deployers,  and  all 
other  panel  one  members  are  presented  in  online  Table  DS1. 
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Notably,  single  and  multiple  deployers  were  proportionately  more 
likely  to  be  male,  active  duty,  air  force  members,  deploy  longer, 
experience  more  combat,  screen  positive  for  binge  drinking,  and 
have  lower  stress  compared  with  the  other  panel  one  participants 
(P<0.05). 

The  examination  of  the  unconditional  models  for  the 
unrestricted  (including  those  with  missing  outcome  or  covariate 
data)  and  restricted  (no  missing  data)  populations  indicated  no 
significant  difference  in  the  resulting  trajectory  solutions.  Fit 
statistics  for  one-  to  five-class  unconditional  models  (no 
covariates)  showed  relative  improvements  in  fit  going  from  two 
to  five  classes  (Table  1).  However,  especially  with  large  sample 
sizes  that  can  accommodate  more  fine-grained  distinctions  among 
classes,  fit  statistics  may  suggest  theoretically  unsound  class 
solutions.  The  four-class  model  for  single  deployers  produced  a 
viable,  theoretically  defensible  solution.  The  five-class  model  was 
unparsimonious  and  unviable  because  it  split  one  class  into  two 
parallel  classes,  creating  a  very  small  class  (1.6%),  and  failed  to 
converge  when  covariates  were  included  in  the  model. 
Accordingly,  we  chose  the  four-class  solution  as  optimal.  Nested 
analyses  of  the  four-class  single-deployer  solution  indicated 
improved  fit  when  the  intercept  variance  was  unconstrained 
(random  effect)  and  when  the  following  covariates  were  included: 
age,  marital  status,  education  level,  service  branch,  pay  grade,  PCS, 
combat  experience,  smoking  status,  heavy  drinking  and  binge 
drinking.  Figure  1  shows  that  the  majority  of  single  deployers  were 
assigned  to  a  class  with  consistently  low  symptoms  (low-stable, 
83.1%).  Other  classes  were  characterised  by  low  pre-deployment 
symptoms  with  elevated  and  worsening  post-traumatic  stress 
after  deployment  (worsening-chronic,  6.7%),  moderate  pre¬ 
deployment  post-traumatic  stress  that  decreased  steadily  after 
deployment  (moderate-improving,  8.0%),  and  elevated  post- 
traumatic  stress  at  all  measurement  points  (high-stable,  2.2%). 
Slope  estimates  were  significantly  non-zero  for  the  worsening- 
chronic  post-traumatic  stress  (28.27,  standard  error  (s.e.)  =  2.05, 
P<  0.001)  and  moderate-improving  (  —  8.86,  s.e.  =  1.24, 
P< 0.001).  The  slope  for  the  low-stable  group  was  small,  but  also 
significant  due  to  the  group’s  large  size  and  small  standard  error 
(1.64,  s.e.  =  0.14,  P<0.001).  The  slope  for  the  high-stable  group 
was  non-significant  (  —  5.07,  s.e.  =  5.50,  P=  0.36). 

Among  single  deployers,  those  with  worsening-chronic,  high- 
stable  or  moderate-improving  post-traumatic  stress  had  a  lower 
PCS  score,  indicating  worse  physical  health  than  the  low-stable 
group  (Table  2).  Interestingly,  heavy  drinkers  had  three  times  greater 
odds  of  experiencing  high-stable  than  low-stable  post-traumatic 
stress,  but  also  threefold  odds  of  moderate-improving  post- 
traumatic  stress  symptoms.  The  worsening-chronic  and  high- 
stable  groups  had  more  combat  experiences  than  the  low-stable 
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Fig.  1  Conditional  model  (including  covariates)  of  post- 
traumatic  stress  over  time  among  3393  Millennium  Cohort 
participants  with  a  single  deployment  between  baseline 
(pre-deployment)  and  first  follow-up. 


Post-traumatic  stress  disorder  (PTSD)  assessed  using  the  PTSD  Checklist  -  Civilian 
Version. 


group.  The  worsening-chronic  group  was  more  likely  to  be  in 
the  army  and  have  a  high-school  education  or  less  than  the 
low-stable  group. 

Similarly  to  the  single  deployers,  the  four-class  unconditional 
model  for  multiple  deployers  produced  a  viable  and  parsimonious 
solution,  whereas  the  five-class  model  again  created  a  very  small 
class  (1.3%)  and  again  failed  to  converge  when  covariates  were 
included.  Thus  the  four-class  solution  was  selected  as  optimal. 
Nested  analyses  of  the  four-class  solution  for  multiple  deployers 
indicated  improved  fit  when  both  the  intercept  and  slope  variance 
were  unconstrained  and  when  the  following  covariates  were 
included:  age,  education  level,  service  branch,  pay  grade, 
combat  experience,  current  smoking  and  heavy  drinking.  Figure 
2  shows  that  the  majority  of  multiple  deployers  were  again 
assigned  to  a  low-stable  class  (84.9%).  Other  classes  included 
worsening-chronic  (4.5%),  moderate-improving  (8.5%)  and 
high-improving  (2.2%)  trajectories.  Slope  estimates  were 
significantly  non-zero  for  the  worsening-chronic  post-traumatic 
stress  (31.35,  s.e.=  1.78,  P<0.001)  and  high-improving 

(  —  15.32,  s.e.  =  2.94,  P<  0.001)  classes.  The  slopes  were  relatively 
small  but  also  significant  for  the  low-stable  (1.29,  s.e.  =  0.10, 
P<0.001)  and  moderate-improving  (  —  4.84,  s.e.  =  0.84, 
P<  0.001)  classes. 


Table  1  Fit  indices  for  one-  to  five-class  unconditional  growth  mixture  models  for  single  and  multiple  deployers 

Fit  indices 

1  class 

2  classes 

3  classes 

4  classes 

5  classes 

Single  deployers  (n  =  3393) 

Akaike  information  criterion 

71  637.96 

69472.81 

68532.31 

67784.48 

66735.33 

Bayesian  information  criterion 

71  693.12 

69546.36 

68611.99 

67882.55 

66864.05 

Sample  size-adjusted  Bayesian  information  criterion 

71  664.52 

69508.23 

68570.69 

67831.71 

66797.32 

Entropy 

0.98 

0.98 

0.97 

0.97 

Lo-Mendell-Rubin  likelihood  ratio  test,  P 

<0.001 

<0.05 

0.29 

<0.05 

Bootstrap  likelihood  ratio  test,  p 

<0.001 

<0.001 

<0.001 

<0.001 

Multiple  deployers  (n= 4394) 

Akaike  information  criterion 

89371.29 

86409.24 

84755.68 

83830.92 

83038.46 

Bayesian  information  criterion 

89428.78 

86485.90 

84851.51 

83945.90 

83172.61 

Sample  size-adjusted  Bayesian  information  criterion 

89400.18 

86447.77 

84803.85 

83888.74 

83105.88 

Entropy 

0.98 

0.98 

0.97 

0.97 

Lo-Mendell-Rubin  likelihood  ratio  test,  P 

<0.001 

<0.01 

0.10 

0.74 

Bootstrap  likelihood  ratio  test,  P 

<0.001 

<0.001 

<0.001 

<0.001 
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Table  2  Adjusted  odds  of  class  membership,  single  and  multiple  deployers3 

Single  deployers,  OR  (95%  Cl)  (n  =  3393) 

Multiple  deployers,  OR  (95%  Cl)  (n  =  4394) 

Classes 

Worsening-chronic 

High-stable 

Moderate- 
im  proving 

Worsening-chronic 

High-improving 

Moderate¬ 

improving 

Physical  component  summary6 

0.95  (0.93-0.97) 

0.91  (0.89  -  0.94) 

0.96  (0.94-0.98) 

Education 

High  school  or  less 

More  than  high  school 

Reference 

0.67  (0.45-0.99)c 

Reference 

0.55  (0.29-1.02) 

Reference 

0.71  (0.52-0.97) 

Reference 

0.66  (0.40  -1.08) 

Reference 

0.70  (0.36-1.34) 

Reference 

0.55  (0.41-0.74) 

Combat  experience6 

1.44  (1.31-1.59) 

1.25  (1.02-1.52) 

1.06  (0.97-1.17) 

1.71  (1.53-1.91) 

1.24  (1.05-1.46) 

1.09  (1.01-1.18) 

Service  branch 

Navy/coast  guard,  marine 
corps,  air  force 

Army 

Reference 

2.57  (1.65-3.99) 

Reference 

0.96  (0.46-2.03) 

Reference 

0.92  (0.66-1.29) 

Reference 

1.78  (1.17-2.69) 

Reference 

0.98  (0.55-1.76) 

Reference 

1.08  (0.82-1.43) 

Pay  grade 

Enlisted 

Officer 

Reference 

0.44  (0.26-0.73) 

Reference 

0.61  (0.25-1.51) 

Reference 

0.42  (0.26-0.70) 

Reference 

0.38  (0.22-0.67) 

Reference 

0.43  (0.18-1.00) 

Reference 

0.51  (0.35-0.76) 

Heavy  drinking 

No 

Yes 

Reference 

1.65  (0.98-2.78) 

Reference 

3.14  (1.57-6.28) 

Reference 

3.08  (2.08-4.58) 

Reference 

1.64  (0.98-2.74) 

Reference 

1 .83  (0.96-3.48) 

Reference 

2.03  (1.41-2.94) 

Binge  drinking 

No 

Yes 

Reference 

1.12  (0.77-1.63) 

Reference 

0.89  (0.40-1.96) 

Reference 

1.25  (0.89-1.74) 

Smoking  status 

Current  smoker 

Never  smoker 

Past  smoker 

Reference 

0.66  (0.42-1.05) 

0.51  (0.30-0.87) 

Reference 

0.43  (0.21-0.85) 

0.58  (0.30-1.10) 

Reference 

0.77  (0.53-1.14) 
0.76  (0.50-1.15) 

Current  smoker 

No 

Yes 

Reference 

1.08  (0.69-1.68) 

Reference 

1.97  (1.22-3.17) 

Reference 

1.71  (1.30-2.26) 

Age6 

1.03  (1.00-1.05) 

1.01  (0.97-1.06) 

0.98  (0.96-1.01) 

1.00  (0.97-1.03) 

0.95  (0.91-0.98) 

1.00  (0.98-1.01) 

Marital  status 

Divorced 

Never  married 

Married 

Reference 

1.89  (0.79-4.51) 
1.59  (0.71-3.55) 

Reference 

1.92  (0.56-6.57) 
1.17  (0.40-3.41) 

Reference 

1.05  (0.57-1.91) 
0.73  (0.42-1.27) 

a.  Empty  spaces  in  table  indicate  that  a  given  covariate  did  not  improve  fit  and  was  not  included  in  the  model. 

b.  Measured  continuously;  relative  odds  based  on  one-unit  increments. 

c.  Bolding  indicates  significance  based  on  a  P<0.05  criterion  in  adjusted  analyses. 

Multiple  deployers  in  the  moderate-improving  and 
worsening-chronic  groups  were  more  likely  to  be  enlisted  service 
members  and  to  report  greater  combat  experience  compared  with 
low-stable  participants  (Table  2).  Participants  in  the  moderate¬ 
improving  class  were  more  likely  to  have  a  high-school  education 
or  less,  report  current  smoking  and  screen  positive  for  heavy 
drinking  when  compared  with  the  low-stable  group. 
Worsening-chronic  post-traumatic  stress  participants  had  greater 
odds  of  being  in  the  army  compared  with  the  low-stable  group. 


Discussion 

Main  findings 

Our  trajectory  analyses  of  the  Millennium  Cohort  data  addressed 
many  of  the  key  methodological  limitations  of  previous  studies  on 
deployed  military  personnel  and  revealed  several  noteworthy 
findings.  First,  across  both  the  single-  and  multiple-deployer 
samples,  the  proportion  exhibiting  a  stable  trajectory  of  low 
post-traumatic  stress  pre-  to  post-deployment,  or  resilience,  was 
exceptionally  prevalent  (83.1%  and  84.9%  respectively).  High 
proportions  of  resilience  have  been  reported  across  a  range  of 
trauma  and  disaster  studies.9,10,14  The  fact  that  the  proportions 
were  even  higher  in  the  current  investigation  was  not  altogether 
surprising.  The  exemplary  resilience  in  these  studies  might  be 
expected,  given  the  level  of  preparedness  and  training  common 
to  US  military  personnel.  Prior  experience  and  training  have 
generally  been  associated  with  favourable  outcomes  in 
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Fig.  2  Conditional  model  (including  covariates)  of  post- 
traumatic  stress  over  time  among  4394  Millennium  Cohort 
participants  with  multiple  deployments. 


Post-traumatic  stress  disorder  (PTSD)  assessed  using  the  PTSD  Checklist  -  Civilian 
Version. 
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populations  at  risk  for  trauma  exposure,27  presumably  because 
they  increase  perceptions  of  control  and  self-efficacy,  help  shape 
constructive  appraisals  and  facilitate  regulation  of  physiological 
and  emotional  reactions.28  Similar  levels  of  resilience  were 
reported  in  the  prospective  study  of  United  Nations  peacekeepers.17 

Extending  previous  population-based  military  studies,  the 
proportion  exhibiting  worsening-chronic  post-traumatic  stress 
was  relatively  modest  (6.7%  of  single  deployers  and  4.5%  of 
multiple  deployers).  Several  factors  highlighted  in  recent 
reviews29,30  help  explain  these  findings.  First,  many  studies  have 
assessed  symptoms  based  on  a  screening  tool  using  four  items 
and  in  proximity  to  deployment,  whereas  these  analyses  are  based 
on  a  17-item  assessment  with  an  average  of  2.7  years  between 
questionnaire  responses.  Second,  the  Millennium  Cohort  is  a 
population-based  sample  including  all  service  branches  as  well 
as  reserve  and  national  guard.  Studies  of  Department  of  Veterans 
Affairs  populations  that  report  higher  levels  of  PTSD  only  account 
for  a  minority  of  veterans  who  self-select  to  receive  treatment.  In 
addition,  elevated  levels  of  PTSD  symptoms  have  been  associated 
with  seeking  medical  compensation,  whereas  the  endorsement  of 
PTSD  symptoms  in  the  Millennium  Cohort  questionnaire  was 
confidential  and  not  tied  to  either  monetary  compensation  or 
treatment.2'31  Finally,  many  studies1,32,33  have  used  retrospective 
methods  to  assess  the  proportion  of  PTSD  symptoms,  which 
may  create  artificial  associations  between  PTSD  and  deployment 
due  to  either  recall  or  selection  bias.  The  prospective  design  of  this 
study  assessed  PTSD  symptoms  at  time  points  before  and  after 
deployment. 

The  final  models  exhibited  the  same  number  and  similar  types 
of  trajectories  for  single  and  multiple  deployers.  The  most  elevated 
post-traumatic  stress  scores  pre-deployment  represented  only 
2.2%  of  the  study  population  in  both  single  and  multiple 
deployers.  However,  post-traumatic  stress  scores  in  multiple 
deployers  showed  greater  improvement  over  time  compared  with 
single  deployers.  In  addition,  fewer  service  members  experienced 
elevated  post-traumatic  stress,  in  general,  among  the  multiple 
deployers.  These  findings  may  indicate  that  at  least  some  service 
members  who  deploy  multiple  times  may  be  able  to  participate 
in  repeated  deployments  because  they  are  healthier  or  perhaps 
they  view  themselves  as  physically  and  mentally  fit  for 
deployment. 

Another  interesting  finding  was  that  small  but  meaningful 
numbers  of  service  members  who  experienced  moderate  or 
elevated  pre-deployment  post-traumatic  stress  showed  some 
improvement  in  post-traumatic  stress  over  time.  This  improve¬ 
ment  pattern  had  been  observed  previously  in  the  United  Nations 
peacekeeper  study,17  as  well  as  in  prospective  bereavement  studies 
involving  end-of-life  caregivers.10,34  It  is  important  to  note, 
however,  that  in  our  study  these  participants,  particularly  the 
high-improving  group,  were  not  entirely  free  of  post-traumatic 
stress  at  follow-up. 

Given  that  the  vast  majority  of  our  study  population  exhibited 
a  resilient  pattern  of  low-stable  post-traumatic  stress  scores,  we 
were  especially  interested  in  predictors  of  this  group.  Compared 
with  those  experiencing  low-stable  post-traumatic  stress,  single 
deployers  with  all  other  stress  patterns  reported  worse  physical 
health  and,  with  the  exception  of  the  moderate-improving  class, 
experienced  more  combat.  The  latter  is  not  surprising  because 
combat  experience  has  consistently  been  demonstrated  as  one  of 
the  strongest  predictors  of  PTSD  following  deployment.1,2,32,35,36 

Among  single  deployers,  the  worsening-chronic  post-traumatic 
stress  class  with  elevated  post-traumatic  stress  post-deployment 
was  also  more  likely  to  be  enlisted  and  in  the  army  compared  with 
the  low-stable  class.  This  is  not  unexpected  since  ground  infantry, 
largely  composed  of  army-enlisted  personnel,  would  be  most  likely 


to  experience  combat.  Among  single  deployers,  service  members 
exhibiting  worsening-chronic  post-traumatic  stress  were  also 
more  likely  to  be  high-school  or  less  educated,  which  is  consistent 
with  lower  pay  grade.  Predictably,  the  worsening-chronic  group 
also  experienced  the  greatest  level  of  combat. 

The  most  prominent  predictors  of  post-traumatic  stress 
improvement  in  both  single  and  multiple  deployers  included 
pre-deployment  heavy  drinking  and  enlisted  status  compared  with 
the  low-stable  groups.  In  multiple  deployers,  current  smoking 
predicted  improvement  and  in  both  single  and  multiple  deployers, 
the  moderate-improving  classes  had  significantly  lower  education 
than  the  low-stable  groups.  Although  still  reporting  post- 
traumatic  stress  symptoms,  it  is  possible  that  these  improving 
groups  were  using  alcohol  or  smoking  as  maladaptive  coping 
mechanisms.  Individuals  reporting  a  pattern  of  pre-deployment 
heavy  drinking  may  be  more  likely  to  use  alcohol  to  self-medicate 
for  post-traumatic  stress  symptoms.  A  recent  study  among 
adolescents  exposed  to  traumatic  events  suggested  that  post- 
traumatic  stress  symptoms  are  associated  with  coping-related 
drinking  motives.37  Additional  factors  that  might  foster  improve¬ 
ments  in  post-traumatic  stress  over  time  are  unit  cohesion  and 
social  support.  Supportive  resources  may  mitigate  the  chronic 
effects  of  combat-related  trauma,  particularly  among  enlisted 
service  members,  because  they  provide  opportunities  to 
participate  in  community-based,  healthcare  provider  or  military 
command-directed  reintegration  programmes. 

Membership  in  the  high-stable  (single  deployers)  and  high- 
improving  (multiple  deployers)  classes  was  predicted  by  similar 
variables,  such  as  being  a  current  smoker.  In  addition,  both 
reported  more  combat  experience  than  their  moderate-improving 
and  low-stable  post-traumatic  stress  counterparts.  Interestingly, 
heavy  drinking  was  not  a  predictor  of  worsening  post-traumatic 
stress  after  deployment  but  rather  associated  with  moderate  and 
elevated  pre-deployment  post-traumatic  stress. 

Strengths  and  limitations 

Using  Millennium  Cohort  data  for  this  study  had  advantages.  It 
provided  a  prospective  assessment  for  both  exposures  and 
outcomes  of  interest,  with  an  average  of  2.7  years  between 
assessments;  it  included  members  from  all  services  and  personnel 
from  the  national  guard  and  reserves;  sample  sizes  were 
sufficiently  large  to  obtain  statistically  meaningful  results;  and 
the  data  allowed  for  adjustment  of  potential  confounders  that  are 
otherwise  difficult  to  obtain  from  sources  other  than  self-report. 

There  were  also  a  number  of  potential  limitations  to  the  study. 
As  there  was  variability  in  the  time  from  post-traumatic  stress 
assessments  to  when  participants  started  and  ended 
deployment(s),  sensitivity  analyses  were  performed  to  test  the 
possible  effects  of  time  from  baseline  (pre-deployment)  to  first 
deployment  and  time  from  last  deployment  to  second  follow-up 
on  the  trajectory  assignments  in  both  the  single-  and  multiple- 
deployer  populations.  When  entered  separately  into  the  final 
models  as  continuous  variables,  these  variables  did  not  increase 
fit  for  either  the  single-  or  multiple-deployer  populations.  We 
concluded,  therefore,  that  the  variability  in  time  did  not  influence 
the  trajectory  solutions  of  the  models. 

In  evaluating  the  data,  we  note  that  although  the  initial 
response  rate  to  invitation  in  the  Millennium  Cohort  Study  was 
36%,  multiple  studies  of  possible  reporting  and  selection  biases 
in  baseline  cohort  data  suggest  reliable  reporting,  minimal 
response  bias  and  a  representative  sample  of  military  personnel, 
including  deployers.2,3,7,38-44  Several  participants  had  missing 
data;  however,  trajectory  solutions  with  and  without  missing  data 
did  not  differ.  It  is  possible  that  the  use  of  self-reported  symptom 
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data  to  identify  individuals  with  chronic  post-traumatic  stress 
might  be  less  precise  than  medical  diagnosis;  however,  previous 
studies  using  questionnaire  instruments  for  PTSD  have  shown 
good  agreement  with  medical  diagnosis.40  Moreover,  in  some  cases 
military  personnel  avoid  seeking  healthcare  for  mental  disorders, 
which  would  likely  lead  to  underreporting  of  post-traumatic  stress 
using  medical  diagnoses.32 

Another  potential  limitation  was  that  because  we  measured 
adjustment  solely  in  terms  of  post-traumatic  stress,  it  is  possible 
that  our  analyses  failed  to  capture  other  forms  of  dysfunction. 
Also,  in  some  cases,  LGMM  analyses  may  identify  non-normality 
rather  than  true  mixture  distributions.  We  consider  this  possibility 
unlikely  given  the  large  representative  sample,  the  high  degree  of 
accuracy  in  classifying  participants  into  trajectories  and  the 
coherent  pattern  of  covariate  prediction  for  each  of  the  trajectory 
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patterns. 

implications 

These  data  offer  the  first  prospective  view  of  post-traumatic  stress 
in  a  large  military  population,  with  measurement  of  post- 
traumatic  stress  prior  to  deployment,  and  indicate  that  the  vast 
majority  of  service  members  remained  without  significant  post- 
traumatic  stress  symptoms  prior  to  and  following  deployments. 
The  next  highest  proportion  of  participants  showed  an  improving 
trend  despite  having  one  or  more  deployments.  Only  a  small 
minority  of  the  participants  exhibited  persistent  or  worsening 
symptoms.  The  most  prominent  predictors  of  these  groups  with 
worsening  or  persistently  high  levels  of  post-traumatic  stress  are, 
for  the  most  part,  difficult  to  change  (for  example  enlisted  army 
service  members  who  experience  combat).  Some  programmes 
have  been  created  and  implemented  to  develop  psychological 
resiliency  (such  as  battlemind  training)  among  these  high-risk 
groups,  however,  perhaps  more  is  needed  in  terms  of  preventive 
interventions  for  PTSD  in  these  specific  groups  of  deployers. 
Smoking,  a  modifiable  behaviour,  also  predicted  membership  in 
many  of  the  groups  with  moderate  or  high  levels  of  post-traumatic 
stress.  This  finding  further  stresses  the  importance  of  programmes 
aimed  at  smoking  cessation  and  the  prevention  of  its  initiation  or 
resumption  among  military  personnel,  particularly  those  who 
experience  higher  levels  of  PTSD.  We  cannot  easily  revisit  previous 
conflicts  to  measure  their  psychological  cost.  However,  the  current 
findings  provide  a  rare  glimpse  at  the  natural  history  of 
post-traumatic  stress  symptoms  in  relation  to  recent  military 
deployments  and  thus  may  also  help  illuminate  the  mental  health 
burden  of  both  past  and  future  conflicts. 
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Data  supplement 


Table  DS1  Pre-deployment  (baseline)  characteristics  of  Millennium  Cohort  Participants3 

Pre-deployment  characteristics 

Panel  lb'c 
(n  =  28  433) 

Single  deployersd 
(fl  =  3393) 

Multiple  deployers 
(n  =  4394) 

Gender,  n  (%) 

Male 

20765  (73.0) 

2696  (79.5) 

3759  (85.6) 

Female 

7668  (27.0) 

697  (20.5) 

635  (14.4) 

Age,  years:6  Mean  (s.d.) 

36.6  (8.9) 

34.9  (8.5) 

33.7  (7.8) 

Education,  n  (%) 

High  school  or  less 

1 1  561  (40.7) 

1435  (42.3) 

1731  (39.4) 

More  than  high  school 

16872  (59.3) 

1958  (57.7) 

2663  (60.6) 

Ethnicity,  n  (%) 

White,  non-Hispanic 

20765  (73.0) 

2473  (72.9) 

3171  (72.2) 

Black,  non-Hispanic 

3045  (10.7) 

312  (9.2) 

375  (8.5) 

Other 

4623  (16.3) 

608  (17.9) 

848  (19.3) 

Marital  status,  n  (%) 

Never  married 

6902  (24.3) 

911  (26.9) 

1155  (26.3) 

Married 

19464  (68.5) 

2258  (66.6) 

2965  (67.5) 

Divorced,  separated  or  widowed 

2067  (7.3) 

224  (6.6) 

274  (6.2) 

Service  component,  n  (%) 

Reserve/national  guard 

12924  (45.5) 

1367  (40.3) 

1614  (36.7) 

Active  duty 

15509  (54.6) 

2026  (59.7) 

2780  (63.3) 

Service  branch 

Army 

13293  (46.8) 

1543  (45.5) 

1873  (42.6) 

Navy/coast  guard 

5645  (19.9) 

507  (14.9) 

494  (11.2) 

Marines  corps 

1090  (3.8) 

156  (4.6) 

234  (5.3) 

Air  force 

8405  (29.6) 

1187  (35.0) 

1793  (40.8) 

Pay  grade,  n  (%) 

Enlisted 

19902  (70.0) 

2503  (73.8) 

3036  (69.1) 

Officer 

8531  (30.0) 

890  (26.2) 

1358  (30.9) 

Occupation,  n  (%) 

Combat  specialist 

5793  (20.4) 

705  (20.8) 

1260  (28.7) 

Healthcare 

3510  (12.3) 

295  (8.7) 

287  (6.5) 

Functional  support/service  and  supply 

8169  (28.7) 

913  (26.9) 

915  (20.8) 

Other  occupations 

10961  (38.6) 

1480  (43.6) 

1932  (44.0) 

Cumulative  days  of  deployment, e'f  mean  (s.d.) 

30.8  (106.2) 

180.0  (117.5) 

366.3  (230.9) 

Combat  experiences/’2  mean  (s.d.) 

0.4  (0.9) 

1.3  (1.6) 

1.7  (1.8) 

Injured  during  deployment/  n  (%) 

No 

28419  (100.0) 

3374  (99.4) 

4354  (99.4) 

Yes 

14  (0.1) 

19  (0.6) 

40  (0.6) 

Smoking  status,  n  (%) 

Never  smoker 

16760  (59.0) 

1972  (58.1) 

2687  (61.2) 

Past  smoker 

7356  (25.9) 

856  (25.2) 

1003  (22.8) 

Current  smoker 

4317  (15.2) 

565  (16.7) 

704  (16.0) 

Binge  drinker,  n  (%) 

No 

16085  (56.6) 

1653  (48.7) 

1993  (45.4) 

Yes 

12348  (43.4) 

1740  (51.3) 

2401  (54.6) 

Heavy  drinker,  n  (%) 

No 

26305  (92.5) 

3099  (91.3) 

4049  (92.2) 

Yes 

2128  (7.5) 

294  (8.7) 

345  (7.9) 

Life  stressors, e’h  mean  (s.d.) 

110.1  (82.9) 

103.1  (78.0) 

96.0  (74.8) 

Physical  component  summary, 6,1  mean  (s.d.) 

53.4  (7.5) 

54.2  (6.7) 

54.7  (6.3) 

a.  Some  numbers  may  not  add  up  to  100%  due  to  rounding. 

b.  Panel  1  refers  to  all  remaining  Panel  1  participants  with  complete  data  who  were  not  classified  as  single  or  multiple  deployers. 

c.  All  individual  characteristics  were  significant  (P<0.05)  except  for  education. 

d.  All  individual  characteristics  were  significant  (P<0.05)  except  for  binge  drinking. 

e.  Variable  assessed  as  a  continuous  variable  in  this  study. 

f.  Assessed  for  the  first  time  on  the  first  follow-up  questionnaire  (2004-2006  survey  cycle). 

g.  The  sum  of  the  number  of  affirmative  responses  to  the  five  individual  types  of  combat  experiences. 

h.  Measurement  based  on  the  Holmes  and  Rahe  Social  Readjustment  Rating  Scale. 

i.  Assessed  with  the  Short  Form  36-Item  Health  Survey  for  Veterans. 
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